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Abstract 

Radio tracking of interplanetary probes is an important tool for navigation purposes as well as 
for testing the laws of physics or exploring planetary environments. The addition of an accelerome- 
ter on board a spacecraft provides orbit determination specialists and physicists with an additional 
observable of great interest: it measures the value of the non-gravitational acceleration acting on the 
spacecraft, i.e. the departure of the probe from geodesic motion. 

This technology is now routinely used for geodesy missions in Earth orbits with electrostatic 
accelerometers. The present article proposes a technological evolution which consists in adding a 
subsystem to remove the measurement bias of an electrostatic accelerometer. It aims at enhancing 
the scientific return of interplanetary missions in the Solar System, from the point of view of funda- 
mental physics as well as Solar System physics. The main part of the instrument is an electrostatic 
accelerometer called MicroSTAR, which inherits mature technologies based on Onera's experience 
in the field of accelerometry. This accelerometer is mounted on a rotating stage, called Bias Re- 
jection System, which modulates the non-gravitational acceleration and thus permits to remove the 
measurement bias of the instrument from the signal of interest. 

The article presents the motivations of this study, describes the instrument, called GAP, and 
the measurement principle, and discusses the performance of the instrument as well as integration 
constraints. Within a total mass of 3.1 kg and an average consumption of 3 W, it is possible to reach an 
accuracy of 1 pm s~'^ for the acceleration measured with an integration time of five hours. Combining 
this observable with the radio tracking data, it is therefore possible to compare the motion of the 
spacecraft to theoretical predictions with an accuracy improved by at least three orders of magnitude 
with respect to existing techniques. 

Keywords Accelerometer; Bias rejection; Non-gravitational acceleration; General Relativity; Fun- 
damental physics; Planetary science 

PACS 04.80.Cc; 06.30.Gv; 07.87.+V 

1 Introduction 

Gravitation is one of the four interactions in the standard model of fundamental physics. General Relativ- 
ity, its curre nt theoretical formulation, is in good agreement with most experimental tests of gravitation 
(Will, 20061) . But General Relativity is a classical theory and all attempts to merge it with the quantum 



description of the other fundamental interactions suggest that it is not the final theory for gravitation. 

General Relativity is also challenged by observations at galactic and cosmic scales. The rotation curves 
of galaxies and the relation between redshifts and luminosities of supernovae deviate from the predictions 
of General Relativity. These anomalies are interpreted as revealing the presence of new components in 
the content of the Universe, the so-called "dark matter" and "dark energy" which are thought to constitut e 
respectively 25% and 70% of the energy content of the Universe ( Copeland et al. . 20061 : Friemanl . 12008 ). 



Their nature remains unknown and, despite their prevalence, they have not yet been detected by any 
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other means than gravitational measurements. Given the immense chahenge posed by these large scale 
behaviors, it is necessary to explore every possible explanation including the hyp othesis that General 



Relativity is not a correct description of gravitation at galactic and cosmic scales (|Aguirre et al.l . 12001 



Noiiri fc OdintsovLl200<T ). 



Testing gravitation at the largest scales reachable by man-made instruments is therefore essential to 

bridge the gap between experiments in the Solar System and astrophysical or cosmological observations. 

The most notable experiment was the test of the gravitation law performed by NASA during the extension 

of the Pioneer 10 & 11 missions. This test resulted in what is now known as the "Pioneer anomaly" 
1 I ' ' — — k 

( Anderson et al . 0998, 2002a), one of the few experimental signals deviating from the predictions of 
General Relativity (iLammerzahl et al. . 2008 : Anderson fc Nietol . 20091 ). 

In a context dominated by the questions of dark matter and dark energy, the challenge raised by 
the anoma lous Pioneer signa l s has t o be faced. Effo r ts have been devoted to the reanalysis of Pio- 
neer data (jMarkwardtl . 120021 ; lOlsenl . 120071 : iLevv et all 120091 : iTurvshev fc Tothl . l2009t ) with the aim of 
learning as much as possible on its possible origin, which c an be an experimental artifact as well as 
a hint of considerable i mport ance for fundamental physics ( Revnaud fc Jaekel 2008t Revnaud et al 



20091 : iTurvshev fc Tothl . l201(1 ). In the meantime, theoretical studies have been devoted to determine 
whether or not the anomalous signal could reveal a scale-dependent modification of gravitation law 

while remaining compatibl e with ot her tests. Among the candidates are in particular metric extensions 

ll Jaekel fc Revnaudl . [2005al |bl. l2006al|bl ) as well as muhifield extensions (|Moffatl . l200"5il2006l : lBruneton fc Esposito-Farese . 
20071 ) of General Relativity. 



Naturally, sever a l mission concepts have been proposed ([Anderson e t al., 2002b; Ditt us et al.l . 12005 



Johann et aU 120081; iBertoIami fc Paramod. 120071; IChristophe et all . 12009, : .Wolf et al., . ,200i) to improve 



the experiment performed by Pioneer 10 fc 11 probes. A key idea in most of these proposals is to embark 
accelerometers so as to measure the non-gravitational forces acting on the spacecraft, whatever may be 
the cause of these forces. The addition of this observable not only improves the accuracy and quality of 
the navigation of the spacecraft but also allows understanding the origin of the Pioneer anomalous signals 
by measuring the non-gravitational acceleration of the spacecraft. Indeed, with such an instrument, no 
models are needed to estimate non-gravitational forces, which removes a source of uncertainty in the orbit 
determination process. 

The main subsystem of the instrument is an electrostatic accelerometer called MicroSTAR, which 
inherits mature technologies. The main evolution of the instrument is to mount MicroSTAR on a rotating 
stage, called Bias Rejection System, which modulates the signal of interest, and thus makes it possible to 
remove the measurement bias of the electrostatic accelerometer. After a presentation of the motivations 
for this study, the accelerometer and the measurement principle will be described. The performance of the 
instrument will then be discussed. Finally, the integration constraints, which depend on the spacecraft 
design, will be briefly presented. A precision of 1 pm s~^ for the acceleration measured with an integration 
time of five hours will be demonstrated. Using these measurements along with radio tracking data, the 
instrument enables to enhance by at least three orders of magnitude the accuracy of the comparison 
between the spacecraft gravitational acceleration and theoretical predictions. 

The instrument studied in this article is an important technological upgrade for the future of funda- 
mental physics in space as well as for the exploration of the outer Solar System. It may indeed be fiown as 
a passenger on planetary missions heading toward outer planets. It will improve the science return with 
respect to objectives in fundamental physics - deep space gravitation tests - as well as in Solar Systeni 
physi cs - study of the gravity field and environment of the outer planets or their moons (jMilani et al 



200ll ) . This idea has been included in the Roadmap for Fundamental Physics in Space issued by ESA in 



2010 (|ESAl . l201f<) . Indeed, these two scientific goals are intimately co nnected since t he gra vitation force 



law is related to the ephemeris of planets, moons and minor objects ( Fienga et al. . 2010( ) as well as to 



the origins of the Solar System (jBlanc et al.l . 120051 ). 



2 Overview and motivations 

The trajectory of a spacecraft, which is tracked by the radio link (e.g. Zarrouati . 1987 : Mover . 2008f l. is 
the result of gravitational and non-gravitational forces. The knowledge of the non-gravitational forces 
makes it possible to compute the geodesic motion of the spacecraft from the radio tracking data. This 
geodesic motion can then be compared to predictions of gravitation theories, allowing to test gravitation 
at the Solar System scale. 
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2.1 Orbit reconstruction without accelerometric measurement 

The following analysis is made in a post-Newtonian framework which is sufficient for the purpose of this 
article. The tracking techniques give access to the Doppler acceleration of the spacecraft along the line of 
sight, aDoppier- This Doppler acceleration includes the non-gravitational acceleration, hng, and a term 
called aother- This last term includes the acceleration of the spacecraft due to the curvature of space-time, 
the interaction of the radio waves with the gravitational field, the speed and acceleration of the Earth 
station tracking the spacecraft, and aberration effects. Calling u the unit vector directing the line of 
sight, this leads to the simplified formula: 

O-Doppler = SNG ' U -|- aother (1) 

When aooppier is the only observable, it is necessary to use models to correct for aNc in order to access 
o-other, which Can be compared to theoretical predicti ons. This app roach h as been imp lemented for all 
interplanetary spacecrafts' orbit determinations so far ( Tortora et al . . 2004; Asmar et al.l . [2005'). But the 



orbit determination is then subject to errors due to uncertainties or inaccuracies in the models. Even 
when the average non-gravitational acceleration is properly corrected, temporal fluctuations are poorly 
known. 

2.2 Orbit reconstruction with accelerometric measurement 

The electrostatic accelerometer with bias rejection GAP measures along two orthogonal axis (directed by 
ei and 62) the absolute value of the non-gravitational acceleration acting on the spacecraft hng- When 
u is in the plane generated by ei and 62, one can retrieve the quantity aother by combining the Doppler 
and the accelerometer measurements. 

This provides an additional observable which enhances the orbit reconstruction for several reasons. 
First, it removes parameters to be fitted in the process and consequently improves the quality of the 
orbit determination. Second, it measures the fluctuations of the non-gravitational acceleration, which 
cannot be taken into account by models. Third, it removes the correlations which appear in the orbit 
determination process between the non-gravitational acceleration and the gravitational acceleration, when 
the former are not rn easured. This approach has been used for CHAMP and GRACE geodesy missions 
( Touboul et al. ,[1991). " n the frame of a test of the gravitation law in the Solar System, the presence of an 



accelerometer becomes essential in order to achieve the best possible accuracy on the orbit determination. 
The instrument presented in this article has been designed for fundament al physics purposes, bu t 



it can also be used for orbit reconstruction around planets or during fly-bys (jChristophe et al.l . 120081 ). 
In terms of planetary scientific objectives, the measurement of the non gravitational acceleration would 
improve the reconstruction of the gravity fields and therefore the knowledge of the body interior structure. 
It would also give information on the atmosphere with spatial and temporal resolutions outreaching by 
far those provided without accelerometric measurement. 

2.3 Sources of non-gravitational accelerations 

In this section, different sources of non-gravitational accelerations are discussed for an interplanetary 
spacecraft. Depending on the power requirement on-board and on the target distance from the Sun 
(among other considerations), the power generation can be made using solar panels or Radio-isotope 
Thermal Generator (RTG). In the first case, the predominant effect is the direct solar radiation pressure 
whereas, in the second one, it is the thermal radiation. 

The direct solar radiation pressure corresponds to the action of solar photons on a unit surface. The 
power carried by solar photons by surface unit at one astronomical unit from the Sun, called do, is approx- 
imately equal to P = 1.366 x 10^ W m~^ (Willson fc Mordvinov, 2003). Considerin g a ballistic coeffi cient 
equals to Cs = 0.1 m^ kg^^, which is the order of magnitude for Laplace mission (|Biesbroekl . |2008[ ). the 



acceleration due to radiation pressure is equal to CbP/c — 4.6 x 10~^ m s~^ at one astronomical unit, 
where c is the speed of light. When the spacecraft is at a distance d from the Sun , this value is mu l tiplied 
by the [d^/dY . The fluctuations of this acceleration have been characterized by Frohlich fc LeanI ( 20041 ) 
for frequencies between 10~^ Hz and 10~^ Hz. 

The anisotropic thermal radiation of the spacecraft induces a force which is difficult to analyze and 
predict. One of the reasons is the time variation of the thermo-optics coefficient s of the materials on th e 
spacecraft and the difficulty to forecast them, as illustrated in the analysis by IBertolami et al.l (|2008l ). 
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Toth fc TurvshevI (|2009l ) or lRievers et"all (|2010[ ). The use of an instrument like GAP removes any problem 



concerning the modeling of these forces since they are directly measured. 

When a spacecraft orbits or fly-bys a planet, the atmosphere induces a force on the spacecraft which 
depends on the relative speed, the density and the ballistic coefficient. Models of atmospheric density 
exist but their accuracy has not been proved and they do not take into account any time variability. 
Therefore, computation of drag forces is plagued by large uncertainties. 

The interplanetary medium also induces a drag force on the s pacecraft during it s cruise. This effect 



has been evaluated for the Pioneer probes and found to be small (jNieto et al.l . 120051 ) 



Finally the magnetic field acts on the spacecraft because of its charge via the Lorentz force. When 
close to a planet, the spacecraft will also experience the radiation pressure due to the planet's albedo. 
Other forces which are not listed here may also play a role (e.g. solar wind, maneuver, gas leakage). 



2.4 Electrostatic measurement 

The electrostatic accelerometer, MicroSTAR, measures a combination of the non-gravitational accelera- 
tions presented above and other terms described later in this article. Along one axis, this external signal, 
which is to be measured, is called s and the actual measurement made by MicroSTAR is called m. The 
measurement chain induces a noise and a bias, respectively called n and &, as well as linear and quadratic 
factors, respectively called 5ki and ^2. The bias corresponds to the low- frequency fluctuations and drifts 
with respect to the duration of the measurement. The relation between these quantities is given by: 

m={l + 6ki)s + k2S^ +b + n (2) 

In the following, 5ki and /c2 will be assumed to be equal to zero without jeopardizing the discussion. The 
purpose of the instrument presented in this article is to remove the bias via a modulation of the external 
signal using a rotation of MicroSTAR with a rotating stage called Bias Rejection System. First, the 
electrostatic accelerometer MicroSTAR is presented, as well as the spectral characterization of its noise. 
Then, the measurement principle and the bias rejection scheme are described. This gives the opportunity 
to give some requirements on the Bias Rejection System needed to achieve the target accuracy. 



3 MicroSTAR: a miniaturized electrostatic accelerometer 



MicroSTAR is an electrostatic accelerometer ( Josselin et al.l . [l999f l based on On era's expertise i n the 

field of acceleromet ry and gravimetry: CHAMP, GRAC E and GOCE missions (iTouboul et al.l, Il999l: 

Marque et al.l . [2008[ ) and the upcoming Microscope mission (|Guiu et al.l . l2007tlHudson et al.l . l2007tlTouboul fc" Rodrigue^ 



20011 ). Ready-to- fly technology is used with original developments aimed at reducing power consumption, 
size and weight. The accelerometer is mounted on a rotating stage, called Bias Rejection System, which 
allows removing the bias and therefore making absolute measurements. An Interface and Control Unit 
(ICU) is in charge of interfacing the instrument with the spacecraft. 

The principle of MicroSTAR is to detect by capacitive measurement any motion of the proof mass 
with respect to the cage and to act, via a feedback control loop, on the proof mass thanks to electrostatic 
forces in order to keep it at the center of the cage. So the proof mass is almost motionless with respect 
to the cage in the control bandwidth of the control loop and at the level of the capacitive measurement 



3.1 Design 

The proof mass is inside a cage composed of six identical plates. This choice has been made to decrease the 
cost and production time of MicroSTAR (Fig. [TJ. All these components are made in the same Ultra-Low 
Expansion (ULE) silica glass manufactured by Corning. The proof mass is a cube whose edge measures 2 
cm and which weights 18 g. It is gold coated with a cutting such that two opposite faces have opposite po- 
larisation. On each ULE pla t es, tw o electrodes are engraved by Ultra-Sonic machining, which is patented 



by Onera (jCampergue et al. ■ Il990l) . The position and attitude of the proof mass with respect to the cage 



(6 degrees of freedom) are controlled by the six pairs of electrodes, such that there is no contact between 
the proof mass and the cage. The electrodes are designed such that degrees of freedom are coupled two 
by two. And the design of the accelerometer is completely symmetric. Consequently, the accelerometer 
cannot be lev itated under 1 g: the condition and performance of the instru ment will have to be tested i n 
a drop tower (|Dittud . [l99ll) . as it is foreseen for Microscope accelerometers (jTouboul fc RodriguesifioOll ). 
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Figure 1: Exploded view of the mechanical core of MicroSTAR. The position and attitude of the proof 
mass is controlled by six pairs of electrodes whose potentials are defined by a control loop. The potentials 
ibVJ, on the proof mass are kept precisely constant by means of reference voltage generators connected 
with two thin gold wires to the proof mass. On each ULE plates, four stops prevent the proof mass from 
touching the electrodes. 




Figure 2: The core of MicroSTAR on the sole-plate with the hermetic housing. The proximity electronics 
is around the housing. 



The core of MicroSTAR is mounted on an hermetic sole-plate, and an hermetic housing ensures 
vacuum and cleanliness. A getter maintains a vacuum smaller than 10^'^ Pa during the whole life of the 
instrument. The proximity electronics is mounted on a Printed Card Board (PCB) around the hermetic 
housing. 

The mass estimate for the accelerometer shown in figure [5] is 1.4 kg. The Bias Rejection System which 
will be added to MicroSTAR shall not exceed 1.1 kg and the Interface and Control Unit 1.0 kg. The total 
mass required to host the instrument on board a spacecraft will therefore be 3.1 kg. 



3.2 Control of the proof mass 

Figure [3] shows in one dimension the principle of the electrostatic accelerometer. The force acting on the 
proof mass due to one pair of electrodes, F, is equal to 



2{e-5xY 



2{e + 5x)' 



(3) 



where eo — 1/{^qc?) is the electric constant ( Mohr et al. . 20081) . S is the surface of one electrode, e is the 
gap between the proof mass and one electrode when the proof mass is at the center of the cage {5x = 0), 
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Figure 3: Schematic of a one-dimensional electrostatic accelerometer. The potentials on the proof mass 
ibVJ, are set via two gold wires and the position of the proof mass 5x is controlled thanks to the potential 
on the electrodes Ve- 



Vp is the potential on one side of the proof mass and Ve is the potential on the electrodes. Assuming that 
Sx <C e, which is a valid assumption when the control loop is running, one obtains 



J-] 
ruA 



OJp^Sx + GeVe + '^SxVe 

Vp 



(4) 



where tua is the mass of the proof mass, Up"^ = 2eo'EVp'^ / {niAe^) is the electrostatic stiffness and Ge = 
—2eoT,Vp/{mAe'^) is the electrostatic gain of the accelerometer. For MicroSTAR, the total force acting 
along one direction is the sum of the forces of the two pairs of electrodes perpendicular to this direction; 
and the torque acting along one direction is proportional to the difference of the forces of the two pairs 
of electrodes. 

Note that, contrary to an harmonic oscillator, this force tends to bring the proof mass out of equi- 
librium because of the positive sign before Sx. Therefore, a control loop is required to maintain the 
proof mass at the center of the cage. When 5x — 0, the force acting on the proof mass is proportional 
to the potential of the electrodes Ve- Therefore, the output of the accelerometer is the quantity GeK- 
The coefficient Ge depends on the value of V^. To keep it perfectly constant, the proof mass is linked 
to a reference voltage source via two gold wires, such that two opposite faces of the proof mass have an 
opposite potential, as in figure |31 

This set up is different from previous accelerometers designed at Onera for which the potential on 
the proof mass was uniformly equal to Vp and the potentials on the electrodes were +Ve and —Ve- This 
choice has been made for MicroSTAR for miniaturization purpose as this new set up allows reducing the 
mass of the electronics. 



Table 1: Parameters of MicroSTAR. 





Size of the proof mass 


20 X 20 X 20 mm^ 


rriA 


Mass of the proof mass 


18 g 


E 


Size of one electrode 


94 mm^ 


e 


Gap between the proof mass and the electrode 


0.3 mm 


Ia 


Inertia of the proof mass 


1200 g.mm^ 



3.3 Electronics 

The electronics is composed of capacitive detectors (one for each couple of electrodes) , a FID (Proportional- 
Integral-Derivative) control and a recombination device. As described above, the outputs of the ac- 
celerometer are the voltage applied on each of the electrodes to control the proof mass. These voltages 
are digitalized through an accurate Analog to Digital Converter of 24 bits and sent to the Interface and 
Control Unit. 

The same electrodes are used for detection and action. In order to separate these two functions, 
the detection is done at 100 kHz by adding to the constant voltage Vp a voltage Vd modulated at this 
detection frequency. The consumption estimate for the electronics is 1.4 W. 
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3.4 Spectral characterization of MicroSTAR's noise 
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Figure 4: Square-root of the Power Spectrum Density of MicroSTAR's noise. The main contributors are 
the gold wires at low frequencies, the capacitive detector at high frequencies and the Analog-to-Digital 
Converter for frequencies between 10~^ and 10~^ Hz. 



In this sec t ion, t he noise of the instrument i s brie fly described. The an alysis relies on m odels 
(jCrassia et all I2OOOI : IChhun et all l2007t [Touboull. l2009(). gro und experiments (jWillemenotl . Il997t ) and 



in-orbit measurements ( Flurv et al.l . 12008 



Marque et al 



2OIOI) . Figure [5] shows, for the parameters given 



in table [Hand for a measurement range equal to 1.8 x 
density of the noise and the different contributions: 



10 ms , the square- root of the power spectrum 



The Proportional-Integral-Derivative controller (PID) is a minor contributor, 
performed using operational amplifiers. 



This function is 



The "Recombination" term corresponds to the noise induced by the coupling of two degrees of 
freedom. The controller computes the translation and rotation which must be applied to the proof 
mass. These quantities then have to be translated into forces that must be applied by the electrodes. 
Some noise is generated in this process by the amplifiers which are used to perform this function. 

The Contact -Potent i al-Diff erence (CPD) refers to voltage fluctuations on the electrodes, also called 
patch effect ( Speakd . 1996 ). The chemical potential, which is the sum of the electrical potential and 
the work function, is constant on the surface. The spatial variations of the electrical potential come 
from the spatial variations of the work function, which is due to atomic composition of the materials. 
In addition, there are temporal fluctuations due to evolution of the material characteristics. 

The "parasitic noise" term includes the several contributions: radiometer effect, asymmetric out- 
gassing, piston noise, thermal radiation pressure noise and the two gold wires, which are the main 
contributors. 

The "Detector" term corresponds to the noise generated by the electronics in charge of measuring 
the proof mass position via the measurement of the capacity of each gap. It is the main contributor 
at high frequencies. 



The "Measurement" term corresponds to the noise of the amplifiers used to measure the potential 
on the electrodes. 
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• The Analog-to-Digital Converter (ADC) introduces a quantification noise in the measurement chain 
which is constant for all frequencies. Its value increases with the range of the instrument and 
decreases with the number of bits available. 

An analytic formula of the square root of MicroSTAR noise power spectrum density as a function of 
frequency can be given for simulation purpose: 



/ \ . ( f 



with 

= 5.7 X 10"" m.s-2.Hz"i/2 (6) 



4 Measurement principle 

This section describes how GAP is used to perform bias-free measurements. First, the link between the 
non-gravitational acceleration and the electrostatic force (cf. equation (HJ) is established. The bias of the 
instruments appears in the relation derived and it is explained how to use the Bias Rejection System to 
remove it. 

4.1 Electrostatic measurement of the non-gravitational accelerations 

The link between the force described by equation (|4]) and the non-gravitational acceleration acting on the 
spacecraft is established in this section. Let us assume that Rq (0,e°,ey,e°) is the reference frame of 
the Solar System, which is supposed to be Galilean. i?i (5*, e^, By, e^) is the reference frame attached to 
the satellite and R2 (C, e^, By, e^) is the reference frame attached to MicroSTAR. O is the barycenter of 
the Solar System, C the center of the cage of MicroSTAR, and S the inertia center of the spacecraft. The 
angular velocity of the frame Ri with respect to the frame Rj is called fii/y The Bias Rejection System 
allows MicroSTAR to rotate with respect to the satellite: it will be assumed that R2 is rotating around 
el with respect to Ri and that the rotation is parametrized by the angle 6. This means that $12/1 = ^6^. 

Below is a list of the forces acting on the proof mass and on the spacecraft. S and A represent 
respectively the satellite and the proof mass. A difference is made between the gravitational forces and 
the non-gravitational forces: 

• F^^g: Non-gravitational force acting on the satellite due to the external environment. It includes 
forces discussed in section [231 

• F^j_j.g: Gravitational force acting on the satellite due to the external environment. This force is 
induced by the mass distribution in the Solar System. 

• F^^^: Non-gravitational force acting on the proof mass due to the external environment. 

• F^j_^^: Gravitational force acting on the proof mass due to the external environment. 

• F^^^: Gravitational force acting on the proof mass due to the satellite. In the following, this force 
will be referred to as self-gravity. 

• Fl"^'^: Force acting on the proof mass due to the satellite via the gold wires. 

• Fl'j^^^: Force acting on the proof mass due to the satellite via the six pairs of electrodes (cf. equation 
®)- 

By applying Newton's second law to the satellite and to the proof mass, one obtains the following 
equation: 

(— + —) F1'4a = —FZf^s + Al - A2 - A3 (7) 

\mA msj ms 

where it is assumed that the proof mass remains at the center of the cage due to the control loop - i.e. 
i?i and i?2 are identical - and the satellite is perfectly rigid - i.e. the distance between C and S remains 
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^f^A (10) 



constant, ms and uriA are respectively the masses of the satelhte and of the proof mass. Ai, A2 and A3 
are equal to: 

Al = Jii/o A SC + J7i/o A (J7i/o A SC) (8) 

A2 = f — + — ) Ff_^ + -J-F^,?_^ + (—^%,^A ~ — F«,^,) (9) 
\mA ms J iriA \mA ms ) 

1 1 

V — 

mA ms. 

In the next section, the two terms Ai and A2 will not be taken into account because they are either 
negligible or known with sufficient accuracy to be removed. This hypothesis will be discussed in sectionO 
Moreover, to simplify the discussion, it will be assumed that r2i/o — 0, i.e. the spacecraft is not rotating 
with respect to the reference frame of the Solar System. 

4.2 Bias rejection 

Let's introduce F^^^'^ and F^^^", the components of F|'f^'=^ in R2; and the 

components of F%°^^ in R2; and F^^^ , F^^ and F^^, the components of F^^^g ^o- 

The electronics of the instruments gives the values of the components of F|?fJ^^ in i?2. However, 
there are defaults in the measurement chain (cf. equation ^ with Ski = k2 — 0) and this translates in 
particular in a bias called "electronic bias". If m^, my and m^ are the outputs of MicroSTAR for each 
axis, the following relations hold 

mc = F2^,f + fee + ^C. Ce{a;,y,z} (11) 
where 6^ is the electronic bias and is the measurement noise. Moreover, equation ([T]) yields 

J^i:r = -F|f + M«^cos(0) + <«sin(0)) (12) 

Fi',r = -F^:J' + ^^{^F,'^^^sini9) + F,l^com) (13) 
peiec ^ -Fif + ^^Fl^ (14) 

where fi — mA/{mA + ms) = mA/ms + o{mA/ms) and 6 is the angle controlled by the Bias Rejection 
System. Considering only the plane perpendicular to the rotation axis, the previous equations lead to 

= cos(0) + Fj);^sin(0)) +6,+n, (15) 

my ^ fi (-Fo^i^ sin(6l) + F^y^ cos{0)) +by + ny (16) 

where = — ^1'°''' (^ ^ {^i 2/}) total measurement bias. The non-gravitational acceleration 

acting on the spacecraft is modulated by the angle 6 whereas the bias of MicroSTAR is not. It is therefore 
possible to separate the bias b^ and by from the true acceleration acting on the spacecraft, using the Bias 
Rejection System. 

One possible way to use these two equations is to perform a spectral separation of the acceleration 
and of the bias. Assuming for example that 9 = ujt, with lo a constant angular speed, m^ and my can be 
multiplied by cos(a;t), after being received on ground. For the x axis, this would lead to the demodulated 
quantity 

= ^i^J!? + f (-^J!? cos(2cjt) + Fo'^y^ sin(2wt)) - (6, + Ffj'*) cos(cjf) - cos(c^t) (17) 

The external acceleration which was modulated at the frequency uj/{2tt) in m^ and my (cf. equations 
p3|) and (HH])) is brought back in the continuous by the demodulation. The bias bx + -Ffx'' supposed 
to vary such that the Fourier transform of the term {b^ + -Ff x'') cos{u!t) is null at low frequencies. By 
applying a low-pass filter to rhx and rhy , one recovers an unbiased measurement of F^^ and F^^ plus 
the low frequencies of the modulated noise n^; cos(a;t). Therefore, the noise is also filtered such that only 
the components around the modulated frequency affect the final measurements. This heterodyne method 
is the simplest way of achieving the recovery of the non-gravitational acceleration. It has however the 
disadvantage of requiring the rotating stage to operate continuously. This increases power consumption 
and may lead to a quicker breakdown of the components. It also requires to have a counter-rotating body 
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Figure 5: Modulation signal used to separate the non-gravitational acceleration from the bias of the 
electrostatic accelerometer. This signal is idealized since the flips from one position to the other are not 
instantaneous. 



in order not to transfer angular momentum to the spacecraft. Finally, the rotation may induce additional 
forces on the proof mass which would spoil the signal. 

A more convenient modulation technique consists in flipping MicroSTAR by tt rad, leading to a signal 
with a period called r. This modulation scheme, shown in figure [5J does not have the drawbacks of the 
previous one. In practice however, the rotation between and tt rad is not as sharp as shown on the 
figure. Indeed, the angular acceleration and angular velocity are kept small enough such that there is no 
loss of control of the proof mass due to forces exceeding the maximum range of control of the electrostatic 
actuation. As it was mentioned in the previous paragraph, the data acquired when the accelerometer is 
rotating may be spoiled by additional forces. Therefore, in the case of the signal considered here, the 
data acquired when the accelerometer is rotating are not used. 

Concerning the precision of the measured acceleration, called Sa, it can be approximated by the 
following formula, whose squared value corresponds to the value of the noise spectrum at the modulated 
frequency divided by the integration time. 




where S'„(/) is given by equation ([S]), T is the integration time and r is the period of the modulation 
signal. This gives only an order of magnitude since the harmonics of the modulation signal are not taken 
into account. An exact formula for da can be derived but is beyond the scope of this article. Assuming 
that T = 10 min, that the integration time is 5 h and that the only sources of noise and errors are the ones 
described in section [XH the precision achieved on the acceleration is then equal to 0.8 pm s~^ according 
to equation The exact derivation of the uncertainty for the signal of figure [S] gives 1 pm s~^, with 

a modulation period r = 10 min, a integration time T = 5 h and a rotating time between the position 
rad and tt rad equal to 120 s. 

It is important to notice that the Bias Rejection System allows rejecting the bias of the two axes 
which are perpendicular to the axis of rotation. To remove the bias on the third axis, another rotating 
stage would be needed. The direction along which the absolute measurement of the non-gravitational 
acceleration is not required must therefore be the direction of the rotation axis. In the case of an 
interplanetary probe, this axis should be perpendicular to the orbit plane. 

4.3 Requirements on the Bias Rejection System 

The Bias Rejection System is similar, in its principle, to filter wheel used for space camera or to rotary 
actuator used for solar array drives, optical mechanism drives, deployment mechanisms and antenna 
pointing mechanisms. This section gives the requirements this subsystem has to meet in order to obtain 
the target accuracy of 1 pm s~^. 
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The rotating stage is controlled by the ICU through a closed loop with an encoder for detection of 
the rotation. The accuracy of the pointing knowledge of the rotating stage, called 6* , is crucial to obtain 
the target accuracy called a*. Calling the external acceleration a^xt, the value of 6* is given by 

0* = arcsin ( —) ~ — (19) 

Assuming an external acceleration equal to 10"^ m s^^ and a target accuracy of 1 pm s~^, the requirement 
on the angle is 0* = 10~^ rad. 

Concerning the wobble of the rotating stage, it has to be smaller than a value called cj)* given by: 

/ 2(7* , 

(j)* = J « V20^ (20) 

V CLext 

With the previous numerical values, the requirement on the wobble is 0* = 4.5 x 10^'^ rad. The square 
root in the value of 0* comes from the fact that only the term at the order 0*^ induces an error on 
the acceleration measured by the instrument (the term of order one appears as a bias and is therefore 
removed) . 



5 Accommodation of the instrument 

In the previous section, a precision of 1 pm s^^ on the acceleration measured by GAP was obtained for 
an integration time of 5 h. Equation ([7]) shows however that the force measured Fl?^^ is proportional 
to the non-gravitational forces acting on the spacecraft F^^^ plus a number of terms, Ai, A2 and A3, 
which can induce uncertainties on the measurement. A3, which describes the gold wire action was taken 
into account, and the two other terms, Ai and A2, were supposed to be negligible. This section discusses 
this hypothesis. Ai describes that the spacecraft is not a Galilean reference frame; and the components 
of A2 are approximately constant in Rq which means that it can be considered as an external acceleration 
which is modulated. 

As soon as the instrument is not at the center of inertia of the spacecraft and the spacecraft is rotating 
with respect to the Solar System reference frame, term Ai is not null. It has not been taken into account 
because it was supposed to be known precisely enough. This knowledge requires to measure with an 
accuracy compatible with the target accuracy of the instrument, the attitude of the satellite, its angular 
velocity and the position of the proof mass with respect to the center of mass of the satellite, i.e. the 
vector SC. The attitude of the satellite and its angular velocity are measured with star trackers. The 
knowledge of SC requires to measure the position of the instrument with respect to the satellite prior to 
launch and to precisely monitor the position of the center of inertia of the spacecraft, which may change 
due to the use of propellant. This knowledge, along with the mass of the spacecraft, is also needed to 
estimate the self-gravity F^^^ which appears in the term A2 

NG 

environment is electro-magnetic because the proof mass is protected by a hermetic housing. Depending 
on the spacecraft on board which the instrument will be, this force may not be negligible. Therefore, it may 
be required to measure the magnetic field in order to remove this contribution from the measurements. 

The difference of the gravitational acceleration can be constrained. By introducing the distance 
from the Sun to the inertia center of the spacecraft and Sr the distance between S and A, the following 
upper bound comes: 



The only non-gravitational force, F^^^, possibly acting on the proof mass due to the exterior 



1 „r3 1 



iG fT*^ 

_ "■" ext-^A _ ex 

rriA ms 



< '-^Sr (21) 



where G is the gravitational constant and Mq the mass of the Sun. Given conservative assumptions - 
the spacecraft is at one astronomical unit from the Sun and 6r is equals to 1 meter - the upper bound is 



s 



equal to 7.9 x 10^ m s^ , which is one order of magnitude smaller than 1 pm 

This short description of the "parasitic" terms in the measurement shows that the accommodation of 
the instrument in a spacecraft requires care. Self-gravity may be reduced and better controlled if the 
instrument is placed at the center of mass of the spacecraft. As for the pointing system and the attitude 
control of the spacecraft, it must be compatible with the accuracy target of the instrument. 
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6 Conclusion 

The electrostatic accelerometer with bias rejection GAP, designed to answer the need for large scale tests 
of General Relativity and for planetary observations, relies on a new concept to make measurements of 
the non-gravitational acceleration with no bias. The main idea is to mount the electrostatic accelerometer 
MicroSTAR on a rotating stage, so that the external signal to be measured is modulated and distinguished 
from the bias of the accelerometer after post-treatment. 

The article presented the design of MicroSTAR as well as its spectral characterization. It also gave 
the main specifications on the Bias Rejection System. For a given modulation scheme, it is possible to 
derive the final precision on the acceleration measurement: with an integration time of five hours and a 
modulation period of 10 min, the precision is 1 pm s^^ on the measured acceleration. 

This opens the possibility to improve the accuracy on the comparison between the motion of the 
spacecraft and theoretical predictions by three orders of magnitude with respect to existing techniques. 
This accuracy leads to new applications for large scale tests of General Relativity as well as for improved 
planetary observations. 
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